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The results from the STAR Collaboration on directed flow (ui), elliptic flow («2), and the fourth 
harmonic (114) in the anisotropic azimuthal distribution of particles from Au+Au collisions at 
^snn = 200 GeV are summarized and compared with results from other experiments and the- 
oretical models. Results for identified particles are presented and fit with a Blast Wave model. 
Different anisotropic flow analysis methods are compared and nonfiow effects are extracted from the 
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data. For V2, scaling with the number of constituent quarks and parton coalescence is discussed. 
For W4, scaling with i;| and quark coalescence is discussed. 

PACS numbers: 25.75.Ld 



I. INTRODUCTION 

In heavy-ion collisions at the Relativistic Heavy Ion 
Collider (RHIC), the initial spatially-anisotropic partic- 
ipant zone evolves, via possible novel phases of nuclear 
matter, into the observed final state, consisting of large 
numbers of produced particles with anisotropic momen- 
tum distributions in the transverse plane. Important 
insights into the evolution may be obtained from the 
study of this azimuthal anisotropy, most of which is be- 
lieved to originate at the early stages of the collision pro- 
cess. Unlike at lower beam energies [1], the measured 
anisotropies at RHIC reach the large values predicted by 
hydrodynamic models and conform to the particle mass 
dependence expected from hydrodynamics in the kine- 
matic region where this type of model is expected to be 
applicable, i.e., for transverse momenta below a couple 
of GeV/c [2]. The large observed anisotropy at RHIC 
is argued to be indicative of early local thermal equilib- 
rium, and the particle mass dependence is highly rele- 
vant to interpretations involving a strongly interacting 
Quark Gluon Plasma phase [2-4]. At larger transverse 
momenta, measurements of azimuthal anisotropy are also 
relevant to the observation of jet quenching [5, 6]. Given 
the current debate around these interpretations, we sum- 
marize star's findings to date in the area of azimuthal 
anisotropy, present additional results for identified par- 
ticles, compare in detail the different analysis methods 
and their systematic uncertainties, compare the data to 
various models, and systematize the results with fits to 
the hydrodynamic motivated Blast Wave model. 

The paper is organized into sections on the Experi- 
ment, Methods of Analysis, Results, comparison of anal- 
ysis methods, comparison of results to various models, 
and Conclusions. The Methods Comparisons section is 
rather technical, dealing with systematic errors, nonflow 
effects, and fluctuations. 



II. EXPERIMENT 

The main detectors of the STAR experiment used 
in these analyses are the Time Projection Chamber 
(TPC) [7] and the Forward TPCs (FTPCs) [8]. The Ring 
Imaging Cherenkov detector (RICH) [9] of the STAR- 
RICH collaboration is also used for particle identifica- 
tion. The cuts on the data for most of the TPC analyses 
are described in Table I, except for the upper pt cutoff 
which often goes higher as shown in the graphs. For the 
FTPCs the pseudorapidity acceptance is 2.4 < \rj\ < 4.2, 
only at least 5 hits are required, the distance of closest ap- 
proach of the track to the vertex (dca) is restricted to less 
than 3 cm, and for the vi analysis the vertex z is opened 



TABLE I: Cuts used in the TPC analysis of Au-|-Au collisions 
at y^SNN = 200 GeV. Vertex refers to the event vertex, fit 
points are the space points on a track, and dca is the distance 
of closest approach of the track to the event vertex. 



cut 



value 



Pt 
V 

multiplicity 
vertex z 
vertex x, y 
fit points 



0.15 to 2.0 GeV/c 
-1.3 to 1.3 

> 10 

-25. to 25. cm 
-1.0 to 1.0 cm 

> 15 



fit pts / max. pts > 0.52 
dca < 2.0 cm 

trigger min. bias 



up to ±50 cm. The RICH detector [9] covers jr;! < 0.30 
with a 20° bite in azimuth. The RICH detector separates 
charged mesons from protons -I- anti-protons identified 
track by track. The admixture of baryons in the meson 
sample is always less than 10%. The momenta of the 
particles identified in the RICH come from tracking in 
the TPC. 

The data were collected with a minimum bias trigger 
which required a coincidence from the two Zero Degree 
Calorimeters, with each signal being greater than 1/4 of 
the single neutron peak and arriving within a time win- 
dow centered for the interaction diamond. The central- 
ity definition, which is based on the raw charged particle 
TPC multiplicity with I77I < 0.5, is the same as used pre- 
viously [10]. The centrality bins are specified in Table II. 
The mean charged particle multiplicity given in the Ta- 
ble II is for the cuts in Table I. The estimated values 
in the Table come from a Monte Carlo Glauber model 
calculation [11]. In this calculation, the number of par- 
ticipants is equal to the number of wounded nucleons. 
The estimated errors shown for the calculated quantities 
come from a linear combination of the changes in the 
quantities caused by reasonable variations in the param- 
eters of the model. Minimum bias refers to to 80% most 
central hadronic cross section. Two million events are an- 
alyzed for this paper. For the analysis involving FTPCs 
only 70 thousand events are available. Errors presented 
for the data are statistical. Systematic errors are mainly 
due to the method of analysis, nonflow effects, and fluc- 
tuations; these will be discussed in Sec. V on Methods 
Comparisons. 

Several methods are used to identify particles. The 
energy loss in the gas of the TPC identifies particles at 
low Pt- For this the probability PID method [12, 13] 
is used requiring 95% particle purity unless otherwise 
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TABLE II: Listed for ^snn ~ 200 GeV for each centrality bin are the range of the % most central of the hadronic cross section 
and its mean value, the mean charged particle multiplicity with its standard deviation spread, the estimated mean number 
of wounded nucleons, the estimated mean number of binary collisions, and the estimated mean impact parameter, with the 
uncertainties in these quantities. 



Centrality bin 



6 



8 



% most central 70 - 80 60 - 70 50 - 60 40 - 50 30 - 40 20 - 30 10 - 20 

(centrality) (%) 73.8 64.1 53.9 44.7 35.2 25.4 15.1 

(M)±a 38±11 76±17 134±24 214±32 323±42 468±53 651±64 

{Nwn)±o- 13±4 26±7 46±9 75±11 114±12 165±12 232±10 

{Nbinary)±(T 11±5 28±10 61±17 120±28 216±38 364±51 587±61 
(6) ± a (fm) 
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5-10 0-5 
7.7 2.3 

819±48 961±56 
298±10 352±6 

825±72 1049±72 



13.2±0.6 12.3±0.6 11.3±0.6 10.2±0.5 9.0±0.5 7.6±0.4 5.9±0.3 4.2±0.3 



2.3±0.2 



stated. In the FTPCs the energy loss is not sufficient 
for good particle identification. The RICH detector can 
separate mesons from baryons up to higher pt. Using 
the characteristic kink decay of , one is able to go to 
higher pt- Strange particles up to high pt arc identified 
by their topological decay. 



For the kink analysis of charged kaons, flow parameters 
were measured for particles that decay in flight within a 
fiducial volume in the TPC. The one-prong decay vertex 
( "kink" ) provides topological identification of the particle 
species with good rejection of background [14] . The main 
sources of possible misidcntification arc pion decays, ran- 
dom combinatoric background, and secondary hadronic 
interactions in the TPC gas. The level of background in 
the analyzed sample was estimated to be 5-10% but is pt- 
dependent. Several cuts were applied to the raw signal 
in order to remove most of the background. Pion de- 
cays were removed by applying a momentum-dependent 
decay angle cut, which exploits differences in the decay 
kinematics. Other cuts were also applied to the dE/dx, 
pseudorapidity, and invariant mass of the parent track 
candidate, to the daughter momentum, and to the dis- 
tance of closest approach associated with the two track 
segments at the kink vertex. Finally, there was a quality 
cut to remove candidates with vertices inside the TPC 
sector gaps where spurious kink vertices can arise [14]. 
Currently the kink method can reconstruct charged kaons 
up to Pt ~ 4 GeV/c. The tracking software has difficulty 
resolving a kink vertex when the decay angle is less than 
about 6°. For kaons with pt > 3 GeV/c, the decay angle 
is almost always around 6° or less, so efficiency falls off 
rapidly above 3 GeV/c. The efficiency also suffers from 
the limited fiducial volume; the kaon must decay inside 
a small sub-volume of the TPC in order to provide ade- 
quate track length for both parent and daughter tracks. 



Other strange particles were identified by their decay 
topology [15, 16]. These methods used for the strange 
particle decays have already been described [15, 17]. 



III. METHODS OF ANALYSIS 

Directed and elliptic flow are defined as the first, wi, 
and second, V2, harmonics in the Fourier expansion of 
the particle azimuthal anisotropic distribution with re- 
spect to the reaction plane. The reaction plane contains 
the collision impact parameter. However, normally mea- 
surements are made relative to the observed event plane, 
and are corrected for the resolution of the event plane 
relative to the reaction plane. The event plane angle is 
defined for each harmonic, n, by the angle, of the 
flow vector, Q, whose x and y components are given by 

Qn cos(n*„) = '^[m cos(n0^)] (1) 

(5„sin(n*„) = y^[wt sm{n(/)i)], 

where the (jji are the azimuthal angles of all the particles 
used to deflne the event plane and the weights, Wi, are 
used to optimize the event plane resolution. In this paper 
the weights for the even harmonics have been taken to 
be proportional to pt up to 2 GeV/c and constant above 
that. For the odd harmonics they have been taken to be 
proportional to rj for \r]\ > 1. 

STAR has previously presented results using differ- 
ent methods of analysis. In the standard method [18], 
denoted by particles are correlated with an event 
plane of the same harmonic. Using this method STAR 
has presented results on elliptic flow (1)2) for charged 
hadrons [5, 19], identified particles [12], strange parti- 
cles [15, 20], and multi-strange baryons [17]. In the N- 
particle cumulant method [21], denoted by Vn{N}, TV- 
particle correlations are calculated and nonfiow effects 
subtracted to first order when A'' is greater than 2. Non- 
flow effects which affect w„ are particle correlations which 
arc not correlated with the reaction plane. Two-particle 
cumulants should give essentially the same results and 
errors as the standard method, but multi-particle cumu- 
lants have larger statistical errors. STAR has presented 
four-particle cumulant results [22] for charged hadrons. 
In three-particle mixed harmonic methods relative to the 
second harmonic event plane, denoted by 'i;„{EP2} when 
n ^ 2, the particles of a different harmonic are corre- 
lated with the well-determined second harmonic event 



5 



plane. With mixed harmonics, nonflow effects are greatly 
suppressed. With this method STAR has reported re- 
sults on directed flow (vi) [23-25] and higher harmonics 
(vi) [23, 26] for charged hadrons. 



A. Directed flow methods 

Because directed flow goes to zero at midrapidity by 
symmetry, the first harmonic event plane is poorly de- 
fined in the TPC. A better way to measure vi is to use 
mixed harmonics involving the second harmonic event 
plane; this also suppresses nonflow contributions at the 
same time. One such method is the three particle cumu- 
lant method which has been described [27]. 

We also measure vi using another mixed harmonic 
technique: we determine two first order reaction planes 
^f^P^i and ^-f^P^^ in the FTPCs and the second or- 
der reaction plane ^fj^^ in the TPC. Using the recently 
proposed notation (see [23]) we denote this measurement 
as wi{EPi,EP2}. 



z;i{EPi,EP2} 



(2) 



$FTPC 



2*JPC 



)) 



FTPC, 



^FTPC. 



2^TPC 



Res(*TPC) 



where the (p of the particle is correlated with the 5'^'^^"-' 
in the other subevent, and 



Rcs(vl/JPC^ 



(C0S[2(*2 - *RP)]) 



(3) 



represents the resolution of the second order event plane 
measured in the TPC. This resolution, as usual, is de- 
rived from the square-root of the correlation of TPC 
subevent planes. For the derivation of Eq. (2) see Ap- 
pendix B. 

This new vi method also provides an elegant tool 
to determine the sign of V2- One of the quantities in- 
volved in the above measurement of ?;i{EPi,EP2} (see 
Appendix B, Eq. (B3) and compare to [18], Eq. (18)) is 
approximately proportional to the product of integrated 
values of and V2- Applying factors for weights and 
multiphcities [18] leads to 



(4) 




B. Elliptic flow methods 

The standard method [18] correlates each particle with 
the event plane determined from the full event minus the 
particle of interest. Since the event plane is only an ap- 
proximation to the true reaction plane, one has to correct 
for this smearing by dividing the observed correlation by 
the event plane resolution, which is the correlation of the 
event plane with the reaction plane. The event plane 
resolution is always less than one, and thus dividing by 
it raises the flow values. To make this correction the 
full event is divided up into two subevents (a,b), and the 
square root of the correlation of the subevent planes is 
the subevent plane resolution. The full event plane res- 
olution is then obtained using the equations in Ref. [18] 
which describe the variation of the resolution with mul- 
tiplicity. 

The scalar product method [22] is a simpler variation 
of this method which weights events with the magnitude 
of the flow vector Q: 



(Q»<,,;(??,Pf)) 



(5) 



where Un.i — cos{n(j)i) + i sm(n(f>i) is the unit vector of the 
i*^ particle. If Q„ is replaced by its unit vector, the above 
reduces to the standard method. Taking into account the 
non-flow contribution, the numerator of Eq. (5) can be 
written as [6, 22]: 

cos2(0pj - = Mv2{pt) V2 + {non-flow} (6) 



where (f>p^ is the azimuthal angle of the particle from a 
given pt bin. The first term in the r.h.s. of Eq. (6) 
represents the elliptic flow contribution, where V2{pt) is 
the elliptic flow of particles with a given pt, and V2 is 
the average flow of particles used in the sum; M is the 
multiplicity of particles contributing to the sum, which 
in this paper is performed over particles in the region 
0.15 <pt < 2.0 GeV/c and \rj\ < 1.0. 

The cumulant method has been well described [21, 28] 
and previously used for the analysis of STAR data [22]. 

To reduce the nonflow effects from intra-jet corre- 
lations at high transverse momentum, we also use a 
modified event plane reconstruction algorithm, where all 
subevent particles in a pseudorapidity region of ]A?7] < 
0.5 around the highest pt particle in the event are ex- 
cluded from the event plane determination. With this 
modified event plane method, the full event plane reso- 
lution is 15-20% worse than with the standard method 
due to the smaller number of tracks used for the event 
plane determination. 



where the index d represents the three detectors used in 
the analysis: FTPCi, FTPC2, and TPC. For each cen- 
trality class Md denotes the corresponding multiplicities 
and Wd are the applied weights (77- weighting for and 
px-weighting for ^'2). 



C. Higher harmonic methods 

Since the second harmonic event plane is determined 
so well, one can try to determine the higher even har- 
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monies of the azimuthal anisotropy by correlating par- 
ticles with the second harmonic event plane. However, 
then the event plane resolution is worse because of the 
various possible orientations of the higher harmonies rel- 
ative to the second harmonic event plane. Taking k 
to be the ratio of the higher harmonic number to the 
event plane harmonic number, and using the equations in 
Ref. [18] we obtain the resolutions in Fig. 1 for ffc2{EP2}. 
This method works when the resolution of the standard 
method (fc = 1) is large and therefore those for the higher 
harmonics are not too low. Also, these fc ^ 1 methods 
use mixed harmonics, which involve multiparticle corre- 
lations, greatly reducing the nonflow contributions. 
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FIG. 2: (color online). Directed flow of charged hadrons 
as a function of pseudorapidity. The measurements of 
ui{EPi, EP2} (circles; centrality 20-60%) agree with the pub- 
lished results of ui{3} (stars; centrality 10-70%). 

has also reported [29] vi values using a two-particle corre- 
lation method. While we approximately agree at 77 = 4.0, 
they have finite values at ij = 2.5-3.0, while ours are close 
to zero, as can be seen for ours in Fig. 2. 



FIG. 1: (color online). The event plane resolutions as a func- 
tion of centrality for i;fc2{EP2}. 

The cumulant method with mixed harmonics has also 
been used for V4 [23]. 



IV. RESULTS 

In the following sections we present results for di- 
rected flow, elliptic flow, and the higher harmonics. Some 
of the graphs have model calculations on them which 
will be discussed in Sec. VI. The tables of data for 
this paper are available at http://www.star.bnl.gov/ 
central/publications/ . 



A. Directed flow, i;i{EPi,EP2} 

The STAR TPC has very good capabilities to measure 
elliptic flow at mid-rapidity, while the FTPCs allow one 
to measure directed flow. Figure 2 plots directed flow 
as a function of pseudorapidity, showing that vi appears 
to be close to zero near mid-rapidity. First, the analy- 
sis was done successfully on simulated data containing a 
fixed vi . For real data, using random subevents in the 
two FTPCs to determine l-f^^^ and ■9™'^^ in Eq. (2), 
the results are in agreement with the published measure- 
ments obtained by the three-particle cumulant method 
wi{3} [23, 24], as shown in Fig. 2. Recently, PHOBOS 
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FIG. 3: (color online). The product of vf and V2. The shaded 
band is the mean value of this quantity with its error, aver- 
aged over centralities 20-60%. Since this quantity is positive, 
elliptic flow is measured to be in- 



The sign of V2 determines whether the elliptic flow is 
in-plane or out-of-plane. Although the sign of V2 had 
been determined to be positive from three particle corre- 
lations [23] , the above new method for vi allows another 
method based on the sign of v\ ■ V2- Since is always 
positive, the sign of ■ V2 determines the sign of V2. 

Averaged over centralities 20-60% we measure vf ■ V2 
in Fig. 3 to be (2.38 ± 0.99) • 10"^ This is only a 2.4 
sigma effect and if 10% systematic errors are assumed 
based on Sec. V for both vi and V2 this becomes a 2.2 
sigma effect. Only the mid-centrality bins are averaged 
because in this centrality region the expected nonflow 
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contributions are much smaller than for the more central 
and peripheral bins. Therefore, with these caveats, the 
sign of V2 is confirmed to be positive: in-plane elliptic 
flow. 



B. Elliptic flow, V2 

There have been many elliptic flow results from RHIC. 
STAR has extensive systematics which we will present 
and compare to the other experiments. Many of the 
graphs will contain Blast Wave model fits which will be 
discussed in Sec. VI D in Model Comparisons. We will 
present data separately for the central rapidity region, 
the forward region, and for high pt- 



1. The central region 




0.20.40.60.8 1 1.21.4 1.61.8 
p, (GeV/c) 

FIG. 4: (color online). Charged hadron V2 vs. pt for the 
centrality bins (bottom to top) 5 to 10% and in steps of 10% 
starting at 10, 20, 30, 40, 50, 60, and 70 up to 80%. The solid 
lines are Blast Wave fits. 

The V2(j>t) values for charged hadrons for individual 
centralities are shown in Fig. 4 with Blast Wave fits per- 
formed assuming that all charged hadrons have the mass 
of the pion. The data are well reproduced by the Blast 
Wave parameterization when pt is below 1 GeV/c. Above 
this limit, the contribution of protons in the charged 
hadron sample becomes significant and changes with cen- 
trality, which challenges the pion mass assumption. Fur- 
thermore it has been found that hydrodynamic flow may 
not be applicable above 1 GcV/c, especially for light par- 
ticles, as new phenomena such as hadronization by re- 
combination may become significant [30]. 

Although all the data presented in this paper were col- 
lected using the full magnetic field (0.5 T) of the STAR 
detector, some data were also collected using half the 
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FIG. 5: (color online). V2 vs. pt for charged hadrons from 
0-50% centrality collisions in comparison to data from PHO- 
BOS [29]. The line is a polynomial fit to the STAR data. 
The gray error boxes represent the PHOBOS systematic er- 
rors. The bottom panel shows the ratio of the PHOBOS data 
to the polynomial fit. 




p, (GeV/c) 

FIG. 6: (color online). «2{4} vs. pt for identified particles 
in the 20-60% centrality range. The solid lines are hydrody- 
namic calculations [31]. 

magnetic field. Below 0.5 GeV/c the half field V2 val- 
ues are lower, especially for the more central collisions. 
These are regions where the V2 values are small. Adding 
the absolute value of 0.0025 to the half field V2{pt) data 
brought the two sets of data into approximate agreement 
in this Pt range. This additive value is for both sets of 
data analyzed with a dca cut of 2 cm as is done in this 
paper. The discrepancy gets worse as the upper dca cut 
decreases. The effect is not understood and none of the 
half field data are included in this paper. However, a pos- 
sible explanation is that the half field data have poorer 
two-track resolution and are more sensitive to track merg- 
ing, giving a negative nonflow contribution. If true, there 
could be a possible small residual systematic effect on the 
full flcld data. However, the V2 results are compared to 
PHOBOS data [29] for 0-50% centrality and < 77 < 1.5 
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FIG. 7: (color online). ■U2{2} vs. pt for charged pions, charged kaons, and anti-protons for the centrality bins (bottom to top) 
5 to 10% and in steps of 10% starting at 10, 20, 30, 40, 50, 60, and 70 up to 80%. The solid lines are Blast Wave fits. 



in Fig. 5. The STAR data is for the TPC integrated also 
for 0-50% centrality. The full field data presented here 
agree well with the PHOBOS data. 

Results from four-particle cumulants, W2{4}, arc shown 
in Fig. 6 for particles identified by energy loss in the TPC. 
Also shown are hydrodynamic calculations [31]. The two- 
particle values, W2{2}, for pions, kaons, and anti-protons 
are shown for the individual ccntralities with Blast Wave 
fits in Fig. 7. We use only anti-protons at low pt due to 
contamination of the proton sample from hadronic inter- 
actions in the detector material. 
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FIG. 8: (color online). V2 vs. pt for particles identified in the 
RICH detector from minimum bias collisions. The lines are 
hydrodynamic calculations [31] for pions (upper line), kaons 
(middle line), and protons (lower line). 

Figure 8 shows V2(j>t) for charged mesons and protons 
+ anti-protons identified in the RICH detector. The ex- 
perimental results are compared to hydrodynamic cal- 
culations [31]. In the hydrodynamic picture, the mass 



ordering of V2 (the lighter particles have larger V2 than 
the heavier particles) is predicted to hold at all transverse 
momenta. Up to pt ^ 2 GcV/c, V2 of charged mesons 
is found to be larger than that of the heavier baryons, in 
agreement with hydrodynamic predictions. Above = 2 
GeV/c, the data seem to indicate a reversed trend where 
the protons + anti-protons might have larger V2 values 
than the charged mesons. 
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FIG. 9: (color online). V2 vs. pt for neutral and charged kaons 
for minimum bias collisions. The K5 values are from Ref. [20] . 
The hydrodynamic model line is from Ref. [31]. The insert 
expands the low pt region to make the kaons from dE/dx 
more visible. 

From the kink analysis the results are shown in Fig. 9. 
There were about 0.4 accepted candidate kaons recon- 
structed per event. 

Results are shown in Fig. 10 comparing STAR data 
for K% and A A out to 6 GeV/c with some PHENIX 
data [32], and with hydro calculations [60]. For kaons, we 
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FIG. 10: (color online). V2 vs. pt for strange particles from minimum bias collisions. The STAR Kg and A + A values are from 
Ref. [20]. The PHENIX data are from Ref. [32]. The hydro calculations are from Ref. [60]. 



can now compare V2(j>t) for neutral kaons, charged kaons 
from kinks, and charged kaons from energy loss identifi- 
cation. This is shown in Fig. 9, where the agreement is 
good, but in the insert one can see that the neutral kaons 
tend to be slightly lower than the charged kaons. 




p, (GeV/c) 



FIG. 11; (color online). V2 vs. pt for charged pions, charged 
kaons, and anti-protons from minimum bias collisions in com- 
parison to similar data from PHENIX. The lines are polyno- 
mial fits to the STAR data. The bottom panel shows the ratio 
of the PHENIX data to the polynomial fits. 

We can also compare our results in more detail at lower 
Pt with those from PHENIX [32]. Figure 11 shows V2{pt) 
for charged pions and anti-protons from the energy loss 



analysis requiring 90% purity, and kaons from the kink 
analysis. The PHENIX results are for \7]\ < 0.35, for 0- 
70% centrality, and for protons and anti-protons com- 
bined. In the Pt range where the data overlap, the agree- 
ment is seen to be good. 

It is interesting to see how azimuthal correlations 
evolve from elementary collisions (p+p) through colli- 
sions involving cold nuclear matter (d-|-Au), and then 
on to hot, heavy- ion collisions (Au-|-Au). A convenient 
quantity for such comparisons is the scalar product. In 
the case of only "nonflow" , the scalar product should be 
the same for all three collision systems regardless of their 
system size. This assumes independent collisions and 
that other effects like short range correlations are small. 
Thus, deviations of the scalar product from elementary 
p-|-p collisions result from collective motion and/or ef- 
fects of medium modification. 

Figure 12 shows the scalar product as defined in Eq. 6 
as a function of pt for three different centrality ranges in 
Au-|-Au collisions compared to minimum bias p-l-p col- 
lisions [6] and d-|-Au coUisions. For Au-|-Au colhsions, 
in middle central events we observe a big deviation from 
p-|-p collisions that is due to the presence of elliptic flow, 
while in peripheral events, collisions are essentially like el- 
ementary p-|-p collisions. The azimuthal anisotropy goes 
up to 10 GeV/c but we cannot distinguish whether it is 
from hydro- like flow or from jet quenching. For pt beyond 
5 GeV/c in central collisions, we again find a similarity 
between Au-|-Au collisions and p-l-p collisions, indicating 
the dominance of nonflow effects. The scalar product in 
d4-Au collisions is relatively close to that from p-f p colli- 
sions but there is a finite difference at low pt . This differ- 
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FIG. 12: (color online). Charged hadron azimuthal correla- 
tions vs. pt in Au+Au collisions (squares) as a function of 
centrality (peripheral to central from left to right) compared 
to minimum bias azimuthal correlations in p+p collisions (cir- 
cles) and d+Au collisions (triangles). The Au-|-Au and p-|-p 
data are from Ref. [6]. 
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FIG. 13: (color online). The top panel shows the charged 
hadron scalar product vs. pt for different centrality classes 
in d-f Au collisions, and minimum bias p-|-p collisions. The 
Au+Au result is put there for a reference. The bottom panel 
shows the ratio of scalar product from d+Au collisions to min- 
imum bias p+p collisions for three different centrality classes. 



ence is small if compared to the difference between middle 
central Au+Au collisions and minimum bias p+p colli- 
sions. If we examine the difference by looking in d+Au 
collisions at different event classes that are defined by the 
multiplicity from the Au side (Fig. 13), wc find that the 
scalar product in d+Au increases as a function of multi- 
plicity class, which is contradictory to Au+Au collisions, 
in which the differences rise and fall as a function of cen- 
trality; a typical pattern that is caused by collective flow. 
The trend in d+Au could be explained by the Cronin ef- 
fect, because in high multiplicity events, the Cronin effect 
is expected to produce more collective motion among soft 
particles in order to generate a high pt particle [33] . To 
further test the Cronin effect hypothesis, we studied the 
asymmetry of the scalar product in d+Au collisions in 
Fig. 14. The ratio of scalar product from the Au side 
divided by that from the dcuteron side is greater than 
one at low pt, and decreases to ~0.9 above 2 GcV/c. 
This indicates that there is more collective motion for 
Pt > 2 GeV/c in the deuteron side and pt < 1 GeV/c in 
the Au side, which is again consistent with the Cronin 
effect. Recently, the Cronin effect has been explained by 
final-state recombination [34]. However the influence of 
recombination on azimuthal correlations needs detailed 
study. In addition to spectra, the scalar product results 
open new possibilities for testing these models. 
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FIG. 14: (color online). The ratio of the scalar product 
from Au side (—1.0 < rj < —0.5) to that from deuteron side 
(0.5 < rj < 1.0) vs. Pt from minimum bias collisions. 



2. The forward regions 

Our measurements of elliptic flow V2{r]) for charged 
hadrons at forward pseudorapidities along with those 
from the central region are shown in Fig. 15. The pub- 
lished results [29, 35] obtained by the PHOBOS collab- 
oration showing a bell-shaped curve are confirmed. We 
observe a fall-off by a factor of 1.8 comparing V2{i] = 0) 
with V2{t] ~ 3). While STAR determined the event plane 
near mid-rapidity, PHOBOS did it at forward rapidities, 
which probably accounts for the slightly less fall-off that 
they see. Both measurements were done using the stan- 
dard method. Figure 16 compares our results for V2 ob- 
tained with the method of two-particle cumulants, V2{2}, 
to that for four-particle cumulants, U2{4}. The difference 
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FIG. 15; (color online). Charged hadron V2 vs. 77 at y^SNN = 
200 GeV for STAR minimum bias (stars) and PHOBOS [29] 
mid-central (15-25%) centrality (circles). The open circles 
are PHOBOS data reflected about mid-rapidity. 



6 't 

5 I- 

4 '-- 

3 I- 

2 I- 

1 '-- 



ccnxoocb V2{2} 




-5 



2 3 4 5 

n 



FIG. 16: (color online). Charged hadron U2 vs. r\ for 20-70% 
centrality collisions including the FTPC regions. The open 
circles are W2{2}, the filled squares are U2{4} and the stars 
are t;2{2} for only FTPC particles, not using the main TPC 
particles. 
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FIG. 17: (color online). Charged hadron W2{2} vs. pt for cen- 
trality 20-70% in the FTPC (2.5 < \r]\ < 4.0) regions (stars) 
compared to U2{4} in the TPC + FTPCs (up-triangles) . The 
down-triangles are Monte Carlo results fit to the up-triangles, 
and the solid circles include the FTPC momentum resolution. 



ward rapidities. Results of embedding charged pions (ne- 
glecting protons) in real Au-|-Au events up to 5% of the 
total multiplicity in the FTPCs were used to estimate 
the momentum resolution as a function of rj and pt- At 
77 = 3.0 the momentum resolution goes from 10% at low 
Pt to 35% at Pt = 2.0 GcV/c, but gets about a factor of 
two worse at 77 = 3.5. In Fig. 17 the MC simulation 
V2{pt) including the momentum resolution of the FT- 
PCs seems to explain the observed flattening by smearing 
low Pt particles to higher pt- Thus we can not conclude 
that the shape of the pt dependence of elliptic flow at 
forward rapidities is different from that at mid-rapidity, 
even though the values integrated over pt are consider- 
ably smaller as shown in Fig. 16. 



3. High Pt 



at mid-rapidity will be discussed in Sec. V. The FTPC 
W2{4} values are not quite symmetric about mid-rapidity, 
but not unreasonable considering the statistical errors. 
Within the errors in the FTPC regions, the values from 
the different methods are about the same. 

Figure 17 shows V2{^}{pt) obtained from the four- 
particle cumulant method. Since there are many more 
particles in the main TPC than in the FTPCs, these 
values are mainly at mid-rapidity. W2{4}, which is much 
less sensitive to nonflow effects, is compared to W2{2} at 
forward rapidities, where nonflow may be small. The ob- 
served flattening at pt values around 1 GeV/c for the 
FTPC measurements might be explained by the momen- 
tum resolution of the FTPCs. To quantify the influence 
of the momentum resolution a Monte-Carlo simulation 
of V2{pt) based on the measurements at mid-rapidity was 
done, but the input rj and pt spectra were obtained from 
measurements of the Au-|-Au minimum bias data at for- 



Hadron yields at sufficiently high transverse momen- 
tum in Au-+-Au collisions are believed to contain a sig- 
nificant fraction originating from the fragmentation of 
high energy partons resulting from initial hard scatter- 
ings. Calculations based on perturbative QCD predict 
that high energy partons traversing nuclear matter lose 
energy through induced gluon radiation [36] . Energy loss 
(jet quenching) is expected to depend strongly on the 
color charge density of the created system and the tra- 
versed path length of the propagating parton. Consistent 
with jet quenching calculations, strong suppression of the 
inclusive high-pt hadron production [10, 37] and back- 
to-back high-pt jet-like correlation [38] compared to the 
reference p+p and d-|-Au systems was measured in cen- 
tral Au-|-Au collisions at RHIC. In non-central heavy- ion 
collisions, the geometrical overlap region has an almond 
shape in the transverse plane, with its short axis lying 
in the reaction plane. Partons traversing such a sys- 
tem, on average, experience different path-lengths and 
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therefore different energy loss as a function of their az- 
imuthal angle with respect to the reaction plane. This 
leads to an azimuthal anisotropy in particle production 
at high transverse momenta. Finite values of V2 were 
measured in non-central Au+Au collisions for pt up to 
~ 7-8 GeV/c [5, 6] using the standard reaction plane 
method and two- and four-particle cumulants. The mea- 
surements of azimuthal anisotropics at high transverse 
momenta with the standard reaction plane method and 
two-particle cumulants are influenced by the contribution 
from the inter- and intra-jet correlations. These correla- 
tions, in general, may not be related to the true reaction 
plane orientation and, hence, are a source of nonflow ef- 
fects. A multi-particle cumulant analysis, which has been 
shown to suppress nonflow effects, may give lower V2 val- 
ues because of the opposite sensitivity of i'2{2} and V2{4} 
to the fluctuations of V2 itself described in Sec. VB and 
Ref. [391. 
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FIG. 18: (color online). V2 vs. pt for charged hadrons for 
different centrality bins. The standard reaction plane method 
is shown by open symbols and the modified reaction plane 
method by solid symbols. 
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FIG. 19: (color online). V2 for the standard reaction plane 
method divided by 112 for the modified reaction plane method 
vs. Pt for charged hadrons in different centrality bins. Er- 
ror bars are not shown as the same dataset is used for both 
methods. 
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FIG. 20: (color online). V2 vs. pt for charged hadrons from the 
modified reaction plane method (solid circles). Open circles 
(from Ref. [6], Fig. 2) are the two-particle cumulant results 
after subtracting the correlations measured in p+p collisions. 
Error bars show statistical uncertainties only. 



Figure 18 shows the differential elliptic flow V2 ob- 
tained with the standard and modified reaction plane 
methods 0jS 0j function of pt for different collision cen- 
tralities. The modified event plane method excludes par- 
ticles within IA77I < 0.5 around the highest pt particle. 
For both methods V2 rises linearly up to pt ~ I GeV/c, 
then deviates from a linear rise and saturates for pt > 3 
GeV/c for all centralities. Fig. 9 shows a similar behav- 
ior. Although the statistical errors are large, we observe a 
systematic difference in Fig. 18 for the V2 values obtained 
with the two methods at high transverse momenta. This 
is better illustrated in Fig. 19, where we show the ra- 
tio of V2 obtained with the standard and modified reac- 
tion plane methods. At low transverse momenta {pt < 
2 GeV/c), the V2 values are very similar for both meth- 
ods. At higher transverse momenta, V2 is systematically 



larger for the standard reaction plane method. For more 
peripheral collisions this effect is larger, and it also be- 
gins at lower pt- The modified reaction plane method 
seems to eliminate at least some of the nonflow effects 
at high transverse momenta (up to 15-20% at pt = 5-6 
GeV/c in the most peripheral collisions). The contribu- 
tion of the azimuthal correlations not related to the reac- 
tion plane orientation has been previously studied using 
p+p collisions [6]. In p+p collisions, all correlations are 
considered to be of nonffow origin. In Fig. 12 the az- 
imuthal correlations in mid-central Au+Au collisions are 
very different from those in p+p collisions in both mag- 
nitude and Pt dependence. Figure 20 shows the modified 
reaction plane results on V2{pt) for charged hadrons of 
centrality 20-60%. We find a very good agreement of 
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V2 from the modified reaction plane analysis with the 
two-particle cumulant results after subtracting the cor- 
relations measured in p+p collisions [6] . Neither of these 
modified methods which seem to be necessary at high pt 
give results which differ from the simple standard method 
below Pt of 2 GeV/c, and thus are not used in the other 
analyses of this paper. 



C. Higher harmonics 

1. The central region 
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FIG. 21: (color online). V2 scaled down by a factor of 2, 
and D4{EP2} vs. pt for charged hadrons from minimum bias 
events. Using a fit to the V2 values, the lower solid line is the 
predicted 114 needed to just remove the "peanut" waist (see 
text). 
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FIG. 22: (color online). Graphs of Vn and Vi^lv\. The dashed 
lines are surface shell Blast Wave fits with no p4 or 54 terms 
(See Sec. VI D) to the charged hadron v-2 minimum bias data. 
The resultant ratio U4/U2 is shown as the lower dashed line 
in the ratio graph (b). The solid lines are the fits with the 
addition of p4 and 54. The resultant ratio 114/112 is shown as 
the solid curve in the ratio graph (b). The dotted line in the 
ratio graph (b) at 1.2 represents the average value of the data. 



Our results for charged hadron 774 and we from this 
study have already been published [23, 26], and v^ijit) is 
shown again in Fig. 21. It also was found that scales 
as v\. The value of Vi^jv\ was found to be 1.2, almost 
independent of pt [26] , as can be seen in the ratio graph 
of Fig. 22 (b). 

Kolb [40] pointed out that for large vi the azimuthal 
shape in momentum space described by the Vn Fourier 
expansion is no longer elliptic, but becomes "peanut" 
shaped. Using our high pt plateau experimental values, 
we show this in Fig. 23. Kolb also gives an equation for 
the amount of vi^ needed to just eliminate the peanut 
waist. Figure 21 shows that the experimental values 
considerably exceed this value. 

Fig. 24 shows the W4{EP2}(pt) values for the individual 
centralities with filled elliptic cylinder Blast Wave fits 
assuming all charged hadrons have the mass of a pion. 

Using the probability FID method [12, 13] for charged 
pions and anti-protons, and a topological analysis 
method for Kg and A -f A , we obtain the W4{EP2}(pt) 
and vi(j)t) values shown in Fig. 25. For pions the v\ scal- 
ing ratio is shown in Fig. 26. To make this graph it was 
necessary to combine data points to get reasonable errors 
bars for the ratio because the U4 values are so small. The 
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FIG. 23: (color online). A polar graph of the distribution 
1 + 2v2 cos(2(^) -f 2v^ cos(40) where <\> is the azimuthal angle 
relative to the positive x axis. Plotted are the distributions for 
V2 = 16.5% showing the waist, va = 3.8% having a diamond 
shape, and both coefficients together. 



resulting scaling ratio is consistent with that for charged 
hadrons shown in the Fig. 22 (b) ratio graph. 
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FIG. 24: (color online). U4{EP2} vs. pt for charged hadrons 
for the centrality bins (bottom to top) 5 to 10% and in steps 
of 10% starting at 10, 20, 30, 40, 50, 60, and 70 up to 80%. 
The solid lines are Blast Wave fits. 



2. The forward regions 

In Fig. 27 the fourth harmonic i'4{EP2} shows an av- 
erage value of (0.4 ± 0.1)% in the pseudorapidity cov- 
erage of the TPC (I77I < 1.2). In contrast, its value of 
(0.06 ± 0.07)% in the forward regions is consistent with 
zero, with a 2cr upper limit of 0.2%. Therefore the rel- 
ative fall-off of V4 from 77 = to 77 = 3 appears to be 
stronger than for V2- This behavior is consistent with 
V4 oc V2 scaling. 



3. High pt 

It has been emphasized that W4 has a stronger poten- 
tial than V2 to constrain jet-quenching model calcula- 
tions [40]. Following the same procedure as described 
in Ref. [6], we plot in Fig. 28 the f4{3} from moder- 
ately high Pt . It should be noted that the two most pe- 
ripheral points go up rather than down as they do for 
V2, in apparent violation of V4/V2 scaling at this high 
Pt- We compare the results with the fourth harmonic 
anisotropy generated by energy loss in a static medium 
with a Woods-Saxon density profile, hard sphere (step 
function in density), and the extreme case - hard shell 
limit. The results are shown in Fig. 28. The dashed 
curve corresponds to the hard shell; the upper and lower 
bands corresponds to a parameterization of jet energy 
loss where the absorption coefficient is set to match the 
suppression of the inclusive hadron yields. The lower and 
upper boundaries of the bands around b = 11 fm corre- 
spond to an absorption that gives suppression factors of 
5 and 3, respectively. Note that compared to the case of 
V2 [6], the calculations are less sensitive to the suppres- 



sion factors (narrow bands). These model calculations 
cannot reproduce the correct sign of W4 over the whole 
range of impact parameters, and neither can they repro- 
duce the magnitude of V4. A similar observation was 
made for the magnitude of V2 in this pt range in Ref. [6] . 
In the present case, evidently the absorption of jet parti- 
cles is not the dominant mechanism for producing U4 in 
this Pt range. 



V. METHODS COMPARISONS 

In addition to the standard and scalar product meth- 
ods already described, there are also several subevent 
methods where each particle is correlated with the event 
plane of the other subevent. If the subevents are pro- 
duced randomly, we will call this the random subs 
method. If the particles are sorted according to their 
pseudorapidity, we will call it the eta subs method. In 
these methods, since only half the particles are used for 
the event plane, the statistical errors are approximately 
•\/2 larger, but autocorrelations do not have to be re- 
moved since the particle of interest is not in the other 
subevent. 

Another method involves fitting the distribution of the 
lengths of the flow vectors normalized by the square root 
of the multiplicity [22, 41, 42]: 



Qn = Q„/vM 



vlM 



dP 



1 



2al 



M 2 ^ 



where /q is the modified Bessel function and 
^2 =0.5(1 + 5„). 



(7) 



(8) 



(9) 



Nonflow effects are fit with the parameter gn- The values 
of M are in Table II. 



A. Comparisons 

To make a precise comparison of the various methods 
we have calculated V2 integrated over pt and 77 for the 
main TPC, and plotted it vs. centrality in Fig. 29 (a). To 
make the comparison valid we have used the same events 
and the same cuts, which are shown in Table I. The in- 
tegrated values have not been corrected for the missing 
regions beyond the integration limits given. The system- 
atic error at the lowest pt values (« 0.2 GeV/c) is prob- 
ably larger than at higher pt , but its contribution to the 
integrated Vn values is small because the yield is so low 
there. For constructing the Q vector, linear pt weight- 
ing was used for all methods except the g-distribution 
method, where no weighting was used. From the agree- 
ment of different software implementations of the same 
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FIG. 25: (color online). ii4{EP2} and 112 vs. pt for identified pions, anti-protons, Kg , and A + A for minimum bias collisions. 
The dashed lines are at 1.2 iif. 
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FIG. 26: (color online). The ratio 114/112 vs. pt for identified 
pions. The dashed line is at V4/V2 = 1.2. 



method, wc estimate a relative systematic error (not in- 
cluded) of at least 2% of the V2 values shown. 

The results fall generally into two bands: those for 
two-particle correlations methods, and those for multi- 
particle methods. The difference is due either to the de- 
creased sensitivity of the multi-particle methods to non- 
flow effects, or to their increased sensitivity to fluctuation 
effects [39] . Thus, the "true" flow values must be between 
these two limits. To expand the graph in order to look 
for small differences we also have plotted the ratios to 
the standard method in Fig. 29 (b). It appears that the 
standard method is about 5% lower than the other two- 



FIG. 27: (color online). Comparison of V2 to i)4{EP2} for 
charged hadrons from minimum bias collisions as a function of 
pseudorapidity. The fourth harmonic (squares) is consistent 
with zero at forward pseudorapidities but not at mid-rapidity. 
V2 is shown by circles, scaled by a factor of 1/3 to fit on the 
plot. The larger dataset available for the TPC only (triangles) 
confirms our measurement of ii4{EP2} at mid-rapidities. 



particle correlation methods. We first thought that this 
might be due to nonflow effects affecting the extrapola- 
tion in the standard method from the subevent resolution 
to the full event resolution. However, it also could be due 
to the fact that the standard method uses twice as many 
particles as the subevent methods, and therefore is less 
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FIG. 28: (color online). High pt charged hadron D4{3} inte- 
grated for 3 < Pt < 6 GeV/c vs. impact parameter b, com- 
pared to models of particle absorption; dashed curve is the 
hard shell, higher narrow band is Woods-Saxon, lower wider 
band is hard sphere. The bands have widths for absorption 
to match the observed range of yield suppression. 



B. Nonflow effects and fluctuations 
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FIG. 30: (color online). Charged hadron V2{2}, «2{4}, and 
W2{6} integrated values as a function of centrality. 
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Particle correlations which are not correlated with the 
reaction plane are called nonflow effects when they af- 
fect Vn- Figure 30 shows the two, four, and six-particle 
integral cumulant V2 values using the cuts in Table I. 
The four- and six- particle results agree, showing that 
nonflow effects are eliminated already with four-particle 
correlations. 
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FIG. 31: (color online). The nonflow parameter, g^, as a 
function of centrality. The solid points are from the cumulant 
method. The open circles are from the g-distribution method. 



FIG. 29: (color online). Charged hadron V2 integrated over 
Pt and rj vs. centrality for the various methods described in 
the text. In panel (b) is shown the ratio of «2 to the standard 
method V2. 



sensitive to nonflow effects. But this does not explain 
why the scalar product method fafls in the band with 
the subevents. The values from the eta subevent method 
decrease for peripheral collisions. This could be due to 
decreased nonflow effects for particles separated in pseu- 
dorapidity. 



Nonflow can be calculated by the difference between 
the squares of the two-particle and four-particle cumu- 
lant V2 values, normalized with the number of wounded 
nucleons from Table II. Thus, [22, 27, 65] 



g2 = NwN-ivl{2}~vl{4}), 



(10) 



which is shown in Fig. 31 for ^snn = 200 GeV, 130 
GeV, and from the SPS at 17.2 GeV [65]. The SPS 
values were divided by the multiplicity used and multi- 
plied by NwN, both given in that paper [65]. From the 
q-distribution method of calculating V2, 52 can be ob- 
tained by the increase in the width of the distribution 
from Eq. (9). (It should be pointed out that in these 
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fits, V2 and (72 are somewhat anti-correlated.) For the 
^-distribution method the g2 values were also divided by 
the multiplicity used and multiplied by Nw^- Thus, all 
four results have been renormalized to use the number of 
wounded nucleons. Instead of being independent of cen- 
trality as originally thought, (72 seems to decrease some- 
what for the more peripheral collisions, but appears to 
have the same shape for all the systems. The 17 GeV 
results may be different from the others because (72 could 
vary with the acceptance of the detector. At 200 GeV it is 
possible that 172 from the (/-distribution method is larger 
than from the cumulant method because of real fluctua- 
tions in V2 broadening the g-distribution. Although the 
definition of q in Eq. (7) removes most of the multiplicity 
dependence of Q, Eq. (8) still contains the quantity M, 
and thus is subject to the spread in M in a centrality 
bin. 




10 20 30 40 50 60 70 80 90 100 

% Most Central 



FIG. 32: (color online). Upper panel: The ratio W2{4}/i;2{2} 
for charged hadrons as a function of centrality. The lines 
are a Monte Carlo Glauber model calculation of £2{4}/e2{2}. 
Lower panel: The nonflow parameter, g2, as a function of 
centrality. The lines are a Monte Garlo Glauber model cal- 
culation of NwN ■ {v2/e)'^ ■ (£2(2} — e2{4}). In both panels 
the solid lines assume nucleons, while the dotted lines assume 
quarks. 



Fluctuations of the true V2 can lead to an increase in 
the W2{2} values and an equal decrease in the W2{4} val- 
ues [22]. In Ref. [39] initial spatial eccentricity fluctua- 
tions are calculated in a Monte Carlo Glauber (MCG) 
model and their possible effect on the determination of 
elliptic flow is estimated. To do this they take 



{e{2}r 



2(e2)2 



where the averages are over events and e is the eccen- 
tricity which will be defined in Eq. (15). The physics 
assumption is that V2 oc s. Figure 32 top panel shows 
e{4}/£{2} for the quark and nucleon MCG. As with non- 
flow, this ratio is smaller than unity over the whole cen- 
trality range, with the largest suppression for the nucleon 
MCG. The 130 GeV [22] and 200 GcV data are in be- 
tween the calculated values, and arc closer to the nucleon 
(quark) MCG results for peripheral (central) collisions. 
When the fluctuations are small it can be shown that 
W2{4} ~ i'2{6}, and from Fig. 30 it is clear that the data 
indeed support this. 

Figure 32 bottom panel shows the calculated 172 due to 
eccentricity fluctuations [39] . In contrast to expectations 
from nonflow, which would predict a constant value of 
g2 vs. centrality, the eccentricity fluctuations reproduce 
the observed drop of about a factor 3 vs. centrality as 
observed in the data. 

Thus it appears that either nonflow or fluctuations can 
explain the two bands in Fig. 29. Most probably it is 
some of both. Since nonflow effects and fluctuations raise 
the two-particle correlation values, and fluctuations lower 
the multi-particle correlation values, the truth must lie 
between the lower band and the mean of the two bands. 
At the moment we can only take the difference of the 
bands as an estimate of our systematic error. 



VI. MODEL COMPARISONS 

This section compares the experimental results with 
model calculations. Measurements of event anisotropy, 
especially elliptic flow V2, are sensitive to the early colli- 
sion dynamics [43-46] . Extracting physics from the huge 
set of presented data is done via a variety of methods, 
ranging from transport models which include really quite 
detailed (and diverse) descriptions of the sub-nuclear dy- 
namics, to hydrodynamic models which make simplify- 
ing assumptions (zero mean free path and thermaliza- 
tion) rendering all dynamic details irrelevant and focus- 
ing all physics on the equation of state. We first con- 
sider schematic concepts like coalescence which propose 
an underlying nature of the flowing constituents and al- 
low observable tests of scaling relations implied by those 
concepts. Finally we use a simple Blast Wave parame- 
terization, which tries to see whether a consistent picture 
of all data can be achieved and to identify what arc the 
required driving features (like geometric anisotropy at 
freeze-out, etc). 



A. Coalescence of constituent quarks 

Models of hadron formation by coalescence or re- 
combination of constituent quarks successfully de- 
scribe hadron production in the intermediate pt region 
(1.5 < pt < 5 GeV/c) [20, 30, 47]. These models predict 
that at intermediate pt , V2 will approximately scale with 
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the number of constituent quarks (n) with vijn vs. ptln 
for all hadrons falling on a universal curve. When hadron 
formation is dominated by coalescence, this universal 
curve represents the momentum-space anisotropy of con- 
stituent quarks prior to hadron formation. This simple 
scaling, however, neglects possible higher harmonics and 
possible differences between light and heavy quark flow. 
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FIG. 33: (color online). Top panel: Identified particle vi from 
minimum bias collisions. The vertical axis and horizontal axis 
have been scaled by the number of constituent quarks in). 
Pions are not plotted. A polynomial curve is fit to the data. 
The possible systematic error is indicated by the gray band. 
Bottom panel: The ratio of v-ijn to the fitted curve. 

Figure 33 (top panel) shows vi vs. pt for the identi- 
fied particle data of Fig. 10, where V2 and pt have been 
scaled by the number of constituent quarks (n) . A poly- 
nomial function has been fit to the shown scaled values. 
To investigate the quality of agreement between particle 
species, the data from the top panel are scaled by the fit- 
ted polynomial function and plotted in the bottom panel. 
For pt/n > 0.6 GeV/c, the scaled V2 of , K"^, p+p, 
and A -|- A lie on a universal curve within statistical er- 
rors. The pion points, however, deviate significantly from 
this curve even above 0.6 GeV/c. This deviation may be 
caused by the contribution of pions from resonance de- 
cays [48]. Alternatively, it may reflect the difficulty of a 
constituent-quark-coalescence model to describe the pro- 
duction of pions whose masses are significantly smaller 
than the assumed constituent-quark masses [30] . 

At the end of Sec. VB we estimated that the V2 val- 
ues from two-particle correlations could be systematically 
high by between about 10 to 20%. This was based on 
the integrated values for charged particles and we do not 
know yet how this varies with pt and particle type. How- 



ever, to indicate this estimated systematic error a shaded 
band of 10% is shown in Fig. 33 (top panel). 
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FIG. 34: (color online). Top panels: The V2 of tt , p, 
and A -I- A from three centrality bins (30-70%, 5-30%, and 
0-5% of the collision cross section) scaled by the number of 
constituent quarks (n) vs. pt/n. Polynomial curves are fit to 
the data excluding the pions. Bottom panel: The ratios of 
V2/n to the fitted curves. 



The V2/n of7r±,p, K% 



and A-f A from three centrality 



intervals are shown in the top panels of Fig. 34. The 
and A + A values are from Ref. [20] . In the bottom 
panels, the ratios to the fitted curves are shown. The 
most central data (0-5%) are thought to be affected by 
nonflow correlations (See Sec. V). For the 30-70% and 5- 
30% centrality intervals, the V2 of p. Kg and A -I- A agree 
with constituent-quark-number scaling for the expected 
Pt/n range above 0.6 GeV/c to within 10%. 

Figure 10 showed that the data for the heavier baryons 
seem to cross over the data for the mesons at sufficiently 
high Pt- The data in Fig. 8 are consistent with this. In 
the low Pt region the heavier particles have lower V2 val- 
ues as expected for the mass ordering from hydrodynam- 
ics. In the intermediate pt coalescence plateau region the 
three quark baryons have a larger V2 than the two quark 
mesons. Thus the experimentally observed cross-over is 
thought to be due to a change in the particle production 
mechanism. 

From a simple parton coalescence model one can calcu- 
late [49] the observed f4/f| scaling ratio in terms of the 
same quantity for the quarks. The relationships between 
meson (M) or baryon (B) v^/v^ and quark {q) W4/w| are 



[V4.lv. 
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l/4+(l/2)[V«l];,/2, (11) 



and 



«l/3+(l/3)[V«l];,/3- (12) 
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TABLE III: The ratio v^jvi for all pt and only for pt/n > 0.6 
GeV/c. 





all pt 


Vtin > 0.6 GeV/c 




1.17 ± 0.01 


1.14 ± 0.02 




1.19 ± 0.04 






3.1 ± 1.2 


2.5 ± 1.3 


P 


1.46 ± 0.53 




A + A 0.97 ± 0.18 


0.87 ± 0.22 



These can be rearranged [49] to relate W4/w| for mesons 
and baryons: 

HvlCs - (2/3) Hvl]ll^, + 1/6. (13) 

The observed W4/w| scaling ratios, which appear to be 
fairly independent of pt in Figs. 22 (b), 25 and 26, 
are shown in Table III. Although in Fig. 33, quark- 
number scaling is shown to work within errors at pt /n > 
0.6 GeV/c for all particles except pions, it appears that 
V4/V2 scaling may be applicable over a wider range of pt. 
Charged hadrons are in the Table but should be used with 
care because they represent a complicated superposition 
of baryons and mesons from different values of pt /n where 
the B/M ratio is strongly dependent on centrality and we 
cannot even assume that the values are a good estima- 
tor for mesons. The kaon values are not accurate enough 
to test the above equation. Even though the pions are 
known to deviate from the constituent quark number co- 
alescence predictions, we can calculate with Eq. 13, from 
the charged pions for the wide pt range, that v^/v^ for 
baryons should be 0.96 ± 0.03. This is compatible with 
the values for anti-protons and A-|- A in Table III. Equa- 
tion 13 would be valuable for testing the concept of quark 
coalescence in an equilibrated medium, but the accuracy 
of the data so far do not allow a conclusion. 

If, in addition, one assumes [49, 50] that the scaling 
relation for the partons is 

vl = {vl)\ (14) 

then from Eq. (11) v^/vl = 1/4 + 1/2 = 3/4. For baryons 
this ratio from Eq. (12) is 1/3-1-1/3 = 2/3, which is even 
smaller. But, one can see in Table III that experimentally 
this ratio is close to 1.2 for charged hadrons and pions, 
so that either the parton scaling relation (Eq. (14)) must 
have a proportionality constant of about 2, or the simple 
coalescence model needs improvement. 

B. Transport models 

Most of the transport model analyses were done for 
charged hadrons, but we will only compare some of the 
models with identified hadrons. Microscopic hadronic 
transport calculations under-predict the absolute ampli- 
tude of V2 by a factor of 2 to 3. However, most of the 



observed features, like mass hierarchies in both the lowpt 
region and the meson-baryon order, are seen in hadron 
transport model calculations [51]. The strength of V2 
should be sensitive to the density and interaction fre- 
quency of the constituents. Indeed, when reducing the 
hadron formation time, the V2 values are found to in- 
crease [51]. In addition, the tests with the parton cas- 
cade models AMPT [52] and ZPC [53] give the correct 
mass hierarchy but require a large parton cross section 
in order to mimic the early development of V2 at mid- 
rapidity. In ultra-relativistic nuclear collisions, hadrons 
may not be the right degrees of freedom to describe the 
early dynamics. At large values of pseudorapidity, how- 
ever, the AMPT [54] model seems able to describe the wi, 
W2, and Vi results without the large parton cross sections 
and string melting. At all pseudorapidities, at the later 
stage, when particle density becomes dilute, transport 
effects will become important [55, 56]. 

For D4 the parton cascade model AMPT [50] with 
string melting and a large parton cross section, does cal- 
culate reasonable values. However, the calculated pro- 
portionality constant in Eq. (14) is about 1, while our 
data with a simple coalescence model [30] imply it to be 
about 2. 



C. Hydrodynamic models 

Azimuthal momentum anisotropies in the final state 
are generated by particle re-interactions from azimuthal 
spatial anisotropy in the initial state. In the hydrody- 
namic framework, these re-interactions are modeled by 
assuming zero mean free-path and therefore local ther- 
malization. Hydrodynamic calculations have been suc- 
cessful at reproducing previously published data on V2 
and spectra [57-59]. 

Hydrodynamic calculations have been shown in Figs. 6, 
8, 9, and 10, with reasonable agreement with the V2 and 
W2{4} data up to pt of 1-2 GeV/c. Additional results 
for V2 at low pt from minimum bias collisions, are shown 
in Fig. 35. Results of K% and A -I- A are from Ref. [20]. 
The hydrodynamic calculations [31, 57, 60] are consistent 
with the experimental results considering the systematic 
errors, such as the matching of the centralities are not 
included. Also, as described in Sec. V, the data could 
be 10 to 20% systematically high. To indicate this in 
the plot a band of 10% of the charged pions is shown. 
The characteristic hadron mass ordering of V2 is seen in 
the low Pt region, where at a given pt, the higher the 
hadron mass the lower the value of V2. This supports 
the hypothesis of early development of collectivity and 
possible thermalization in collisions at RHIC [57, 59], 
although the underlying mechanism for the equilibration 
process remains an open issue. 

As seen in Fig. 33 the observed values of V2 saturate 
and the level of the saturation seems dependent on the 
number of constituent quarks (n) in the hadron. The 
saturation value is about 0.07n for pt/n> 1 GeV/c. Hy- 
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FIG. 35: (color online). V2{pt) for charged tt, Kg , p and 
A + A from minimum bias collisions. Hydrodynamic calcu- 
lations [31, 60] are shown as dot-dashed lines. The possible 
systematic error is shown at the bottom. 

drodynamic calculations do not saturate in this pt region. 
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FIG. 36: (color online). Charged n plots on the left and p 
plots on the right for V2 for three centrality bins are shown as 
a function of pt. The data are from centralities 40-50% (open 
triangles), 20-30% (open squares) and 0-5% (open circles). 
The corresponding results of a hydrodynamic calculation are 
shown as dot-dashed lines, solid-lines, and dashed-lines, re- 
spectively. Plots on the top are for a — and plots on the 
bottom are for a — 0.02 fm~^. Here a determines the initial 
velocity kick for the hydrodynamic model calculation [61]. 

Figure 36 shows the centrality dependence of pion and 



anti- proton V2 compared with hydrodynamic results [61]. 
The three centrality bins shown are described in Table II. 
Systematic uncertainties, such as the matching of the 
centralities, are not included. Also, from the Fig. 29 (b) 
ratio graph in Sec. V it can be seen that the 0-5% cen- 
trality data could be 25% high. An important concern 
for the 0-5% centrality bin is the fluctuations. Just av- 
eraging over the spread in impact parameters in this bin 
could lower V2 a factor of two [22, 39]. In the hydrody- 
namic calculation, the decoupling temperature was set to 
100 MeV. In order to fit the pt spectra of (anti-)protons, 
the hydrodynamic evolution was started with an initial 
transverse velocity kick of tanh(a • r), where a is a pa- 
rameter [61]. The results for V2 are shown in Fig. 36. For 
a = 0, Fig. 36 (a) and (b), neither pion nor anti-proton 
results can be fitted . For a = 0.02 fm~^. anti-proton 
(d) V2 can be fitted reasonably well but, for pions (c), 
the model results still miss the data. It appears that 
with the initial velocity, there is too much kick for pions 
at both mid-central and central collisions. Due to their 
light mass, perhaps pions decouple from the system rel- 
atively earlier than protons, as also indicated in the pion 
interferometry results [62] . It seems that for the 40-50% 
centrality data the hydro calculations over-predict the 
data, which is not surprising for peripheral collisions. 

Both Hirano [63] and Heinz and Kolb [64] explain the 
fall-off of V2 at high 77 as being due to incomplete ther- 
malization. The particle density, dN/dy, also falls off in 
the same way, and at high rj is similar to that at mid- 
rapidity at the SPS [65], where the flow values are also 
lower. Possibly, the lower particle density leads to less 
thermalization, and therefore smaller V2 values. 

Hydrodynamic inspired fits have been done for spec- 
tra [66]. Csanad et al. now report results where the 
authors claim that the resulting pt spectra, interferome- 
try parameters, and anisotropy can all be fitted [67]. In 
particular, they have a fall-off of V2 at high rj. But their 
^1(77) has a large wiggle near mid-rapidity which is not 
observed. They further determined the source parame- 
ters and concluded that about 15% of the hadrons are 
emitted directly from the super-heated region. 

So far there have been very few model calculations of 
V4. However, the magnitude and even the sign of U4 are 
more sensitive than V2 to initial conditions in the hy- 
drodynamic calculations [40] . This calculation predicted 
V4/V2 to vary from 0.7 to 0.3 going from low to high pt, 
which is about a factor of two lower than observed in 
the Fig. 22 (b) ratio graph and Table HI. This calcula- 
tion also predicted a strongly negative wg, which is not 
observed [23]. 



D. Blast Wave models 

Blast Wave models parameterize the coordinate and 
momentum freeze-out configuration generated in hydro- 
dynamic calculations. In a self-consistent hydrodynamic 
calculation, this configuration is determined by the Equa- 
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tion of State and freeze-out prescription; in Blast Wave 
calculations, parameters of the distribution may be var- 
ied arbitrarily to fit the data. In this sense, Blast Wave 
is a "toy" model useful mainly to characterize the data 
and determine the magnitude of thermal (random) mo- 
tion, collective motion, geometry, etc. The model also 
provides parameters that can be used to study the evo- 
lution of flow varying the initial conditions, which in this 
paper is achieved by varying centrality. 

The present paper uses two versions of the Blast Wave 
parameterization. In the first one, all particles are emit- 
ted from a surface shell boosted by a constant flow ve- 
locity [12, 26]. In the second one, particles are emitted 
from a filled elliptic cylinder boosted perpendicular to 
the surface of the cylinder, and with a linear transverse 
rapidity profile inside the cylinder [68]. In this paper, 
unless otherwise specified. Blast Wave fits have referred 
to the filled elliptic cylinder version. 
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FIG. 37: (color online). Schematic illustration of an elliptical 
sub-shell of the source. Here, the source is extended in the 
direction out of the event plane {Ry > Rx). Arrows repre- 
sent the direction and magnitude of the flow boost. In this 
example, p2 > 0. From Ref. [68]. 

In recent versions of Blast Wave models, the system is 
assumed boost-invariant in the beam direction. As sug- 
gested in Fig. 37 for the filled elliptic cylinder, the geom- 
etry in the transverse direction is a filled ellipse with the 
major axis aligned with the reaction plane or perpendic- 
ular to it. One may quantify the geometrical anisotropy 
of the system with the eccentricity 
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where the x direction is in the reaction plane. Superim- 
posed on a randomly-directed energy component quan- 
tified by a temperature, T, each geometrical cell of the 
system is boosted "outward" by a velocity (flow) field. 
Here, "outward" indicates the direction normal to the 
surface of the elliptical shell on which the element sits. 
The magnitude of the flow field vanishes (by symmetry) 
at the center of the system and grows linearly with the 



distance from the center, reaching its maximum at the 
transverse edge of the system (here assumed to be a 
sharp, non-diffuse edge). The average value of the fiow 
magnitude is quantified by a parameter pQ. The flow 
magnitude may be larger (or smaller) for sources emit- 
ting in the x- versus the y-direction; the magnitude of 
this boost oscillation with azimuthal angle is quantified 
by the parameters p2 and ^4. In Fig. 37, a larger in- 
plane than out-of-plane boost (corresponding to p2 > 0) 
is suggested by the longer boost angles in-plane. 

Several parameters of the system affect V2- Obviously, 
the larger the magnitude of p2, the larger the momentum- 
space anisotropy. Further, the geometric anisotropy 
plays a role even if the boost strength is identical in all 
directions {p2 = p4 = 0), ii Ry > Rx {Ry < Rx) it is 
clear from Fig. 37 that a greater (lesser) number of el- 
ements boost particles into the reaction plane, resulting 
in anisotropy in azimuthal momentum space. Finally, it 
is clear that the temperature, T, plays a role, since if the 
random energy component is dominant (T larger than 
the rest mass), momentum anisotropics will be reduced. 
An extensive discussion of the interplay between these 
effects may be found in Ref. [68]. 

To summarize, the free parameters of the fits in the 
shell case are T, po, P2, Pi, S2 and S4, where T is the tem- 
perature parameter, the pn are the harmonic coefficients 
of the source element boost in transverse rapidity, and 
the Sn are the harmonic coefficients of the source density 
which boosts into a particular direction. In previous pa- 
rametcrizations [12] where there was no p4, p2 was called 
Pa. In the filled ellipse case the free parameters are T, 
Po I P2 J P4 , Rx J and Ry , where Rx is the in-plane radius of 
the ellipse and Ry is kept constant at a non-zero value. 
In fitting data with a surface shell model po is about 2/3 
as large as for a solid cylinder with a linear profile. The 
eccentricity is approximately equal to 2s2. For an ellipse, 
the parameter S4 is approximately equal to Sj. The ac- 
tual equations used are given in Appendix A. 

First we verified that the hydrodynamic calculations 
reported in Ref. [40] can be successfully fit by the Blast 
Wave model with reasonable parameters: T = 93 MeV, 
Po = 0.91, p2 = 0.080, Pi = 0.0017, e = 0.122. Since the 
hydro had no error bars there is no 'x^ /ndf . While spec- 
tra and V2 are well reproduced up to pt = 1.5-2 GeV/c, 
the Pt dependence of V4 appears quadratic in the Blast 
Wave, while rather linear in the hydrodynamic calcula- 
tion. 

We have seen in Fig. 7 that the Blast Wave parame- 
terization does a good job at simultaneously reproducing 
pion, kaon and anti-proton V2 ■ The fits are performed si- 
multaneously to spectra as well as on V2 and 1)4, in order 
to be over-constrained. Pion, kaon and proton spectra 
(not shown) are well reproduced. Because spectra have 
typically more data points and smaller error bars, both 
T and po can be determined, while p2, Pi, and e are con- 
strained by the ?;„. The total per degree of freedom 
varies for different centralities around an average value of 
56/65; without exhibiting any specific dependences. The 
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average ^ per data point is 14/6 for pions, 7/4 for kaons 
and 17/10 for protons. When fooking at individual data 
sets (e.g pion vi, proton spectra), the is compared to 
the number of data points because the degrees of free- 
dom can only be calculated including all the data points 
as each parameter is constrained by more than one data 
set. Because the error bars are small (less than 5%) 
compared to the spectra error bars (between 5 and 10%) 
the total is dominated by the contributions from the vi 
results. The calculation fits the peculiar negative values 
of the anti-proton V2 in Fig 7 (c) in central collisions with 
•pt below 0.5 GeV/c. This feature is reproduced when p2 
is significant while the thermal velocity is small. In this 
case the flow boost is strong enough that it suppresses the 
low "Pt anti-proton emission in-plane compared to out-of- 
plane [31]. When the eccentricity is sufficiently large this 
phenomenon does not take place. The pion V2{pt) data 
points in Fig 7 (a) are similar in the three most periph- 
eral bins. However, the anti-proton values arc not, and 
thus meaningful fits are still possible. The pt ranges in 
GeV/c used for the Blast Wave fits where the data had 
reasonable error bars were 0.4 to 1.0 for pions, 0.15 to 
0.5 for kaons, and 0.3 to 1.1 for anti-protons. 

The Blast Wave parameters obtained from fitting V2 
and Vi data are shown in Fig. 38. They provide a good 
way to systematize a large amount of experimental data. 
It should be emphasized that other formulations of the 
Blast Wave model would give different fit parameters [69] . 
As the parameters T and po arc constrained mostly by 
spectra, they agree with the values published [70]. p2 
and e are fully constrained by the V2 data. p2 reaches a 
maximum in the centrality region 30-60%. This is easily 
understood recalling that in this centrality region, the ini- 
tial spatial azimuthal anisotropy of the system is large, 
while the initial energy density is still large enough to 
trigger a significant collective expansion. This expansion 
is clearly visible comparing the initial and final eccentric- 
ities. The system spatial deformation is a maximum in 
the region where the azimuthal push quantified by p2 is a 
maximum. Thus, the Blast Wave parameterization pro- 
vides an intuitive self-consistent description of the data. 

For one centrality we show in Fig. 39 the charged 
hadron results from this standard event-plane analysis, 
together with pion results for a standard analysis and a 
two-particle cumulant analysis. As shown in Sec. V, the 
integrated two-particle cumulant U2{2} values are usu- 
ally 5% higher than the standard V2 values. The charged 
hadron values are somewhat smaller than the pion val- 
ues, because of the presence of protons. Even though the 
flow values are fairly close, the p2 fit parameters shown 
in Fig. 38 differ appreciably. This is because the e val- 
ues come out the same and the small differences in the 
V2 values are all forced into the p2 values. It appears 
that the e values are at least half as large as the initial 
eccentricities of the overlap region. 

Both hydrodynamic calculations and Blast Wave fits 
can well reproduce transverse momentum spectra and 
second- harmonic anisotropy (u2). However, as mentioned 
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FIG. 38: (color online). The Blast Wave parameters p2, pi, 
e, po and T plotted vs. centrality. The circles are for pions 
from a two-particle cumulant analysis, the squares for pions 
from a standard event-plane analysis, and the the triangles 
for charged hadrons from a standard analysis. The lines are 
polynomial fits. In the middle panel the initial geometrical 
eccentricity is also plotted as a dashed line. The actual pa- 
rameter values are available at http://www.star.bnl.gov/ 
central/publications/ . 



above, hydrodynamic calculations do not agree with mea- 
sured values of V4. The question, then, is whether Blast 
Wave parameters may be adjusted to simultaneously fit 
V2 and W4, hopefully providing useful feedback to theo- 
rists doing the hydrodynamic calculations. Blast Wave 
fits to V4 are shown in Fig. 24. 

Even with only second-harmonic anisotropics in 
flow strength and spatial geometry, fourth-harmonic 
momentum-space anisotropies (^4) are produced in Blast 
Wave calculations. Thus it is possible that one could gen- 
erate V4 without any fourth harmonic anisotropy pi in 
Eq. (Al) [49]. Using the surface shell Blast Wave model, 
we have fit the V2 data using only p2 and S2, and then 
calculated V4 as shown in Fig. 22 as the dashed lines. The 
calculated V4 values are much too small, indicating that a 
real fourth harmonic term is necessary. Then we allowed 
P4 and S4 to vary as well and obtained the fits shown 
in Fig. 22 as solid lines. The fact that these parameters 
arc significant suggests that the spatial distribution of 
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p, (GeV/c) 

FIG. 39: (color online). For centrality 20-30% we show V2{2} 
for pions (circles), V2 for pions (squares), and V2 for charged 
hadrons (triangles). The solid lines are Blast Wave fits. 



the system initial state has a significant fourth harmonic 
component, which translates into a fourth harmonic flow 
oscillation. 



VII. CONCLUSIONS 

All the presently available STAR data for anisotropic 
flow in Au+Au collisions at ^snn = 200 GeV are pre- 
sented for charged particles and for identifled species. 
Agreement between flow data for STAR and other RHIC 
experiments is good. New evidence confirms our earlier 
finding that elliptic flow is in-plane at RHIC. V2 as a 
function of pseudorapidity is not flat, but conflrmed to 
be bell shaped. A detailed comparison of flow analysis 
methods, shows that either nonflow effects or fluctuations 
can explain the difference between V2 from two-particle 
correlation results and multi-particle correlation results. 
The mass dependence of V2 at low pt follows the pattern 
predicted by hydrodynamic models, but a transition to 
a behavior consistent with quark coalescence at higher 
Pt is observed. For identified particles, V2 scales with 
the number of constituent quarks, n, within errors above 
Pt/n ^ 0.6 GeV/c for charged and neutral kaons, for 
anti-protons, and for A + A hyperons. This supports 
the picture of hadron production via coalescence of con- 
stituent quarks involved in collective anisotropic motion. 
If confirmed it would be a strong argument for the decon- 
finement reached in the system. Only pions deviate from 
this behavior, which partially can be explained by the 
large resonance decay contribution to pion production, 
and by the light pion mass. For the higher flow harmonics 
of order n, Vn scales with V2^'^, consistent with quark coa- 
lescence. However, the ratio V4/V2 is unexpectedly large. 
Some hadronic transport models are a factor of 2-3 lower 



than the data, but others achieve reasonable agreement. 
However, hydrodynamic model calculations provide the 
best predictions for V2 compared with data. The charac- 
teristic collectivity feature - hadron mass dependence in 
the low Pt region - is observed. Hydrodynamic models 
seem to work for minimum bias data but not for central- 
ity selected pion and anti-proton data. The discrepancy 
for the central collision data may be due to nonflow ef- 
fects and fluctuations in the data, and for the peripheral 
collisions from a failure of hydrodynamics. Perhaps, more 
work is needed to improve the hydrodynamic fits, espe- 
cially for the different centrality bins, in order to make 
the case for early thermalization of collisions at RHIC. 
Awaiting further theoretical input or explanation are a 
number of STAR results, such as the large V2 at high 
Pt [6] and the V4 observations, is highly sensitive to 
initial conditions and the equation of state used in hy- 
drodynamic calculations, and therefore a challenge to all 
model descriptions. The data were systematized with 
fits to a Blast Wave model. The Blast Wave framework 
is capable of describing the large volume of experimental 
data up to pt of 1 or 2 GeV/c using a relatively small set 
of fit parameters in each centrality interval, and the fit 
parameters are found to vary smoothly with centrality. 
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APPENDIX A: BLAST WAVE EQUATIONS 

In both the surface shell and the filled elliptic cylin- 
der cases the transverse rapidity parameterization is 
extended to account for a possible fourth harmonic 
anisotropy: 

p{(j)b) = po + P2 cos(20b) 4- p4 cos{4:(j)b), (Al) 

where the flow magnitude and anisotropy arc accounted 
for by the p„ parameters and (jjb is the azimuthal angle 
of the boost source element defined with respect to the 
reaction plane, as shown in Fig. 37. 

The distribution of source elements relative to in 
the case of a surface shell is written including 2"^^ and 
4"^ harmonic azimuthal anisotropy quantified by the S2 
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and S4 parameters respectively: 

fl{(pb) = 1 + 2s2 cos{2(f>b) + 2s4 cos(4(/)6) (A2) 

When S2 is positive more particles are boosted in-plane 
than out-of-plane. In the case of a filled ellipse, the 
boost direction {cj)t) is assumed to be perpendicular to 
the freeze-out surface, which leads to a relationship be- 
tween the space and boost azimuthal angles of the emit- 
ted particles 

tan((/),) = [Ry/Rxf tan((/.fc), (A3) 

with Rx and Ry the in-plane and out-of-plane radii, re- 
spectively. For the w„ analysis R = {R"^ + Ry)^^^ is an 
arbitrary radius, but when interferometry data are also 
fit. the units become significant. The system is bounded 
within an ellipse such as il{r^(f>s) — 6(f((j)s)) with the 
step function and 

= ^{r cosies) /RxY + {rsin{4>s)/RyY. (A4) 
I 



In the filled ellipse case there is no explicit 2"*^ and 4**^ 
harmonic parameterization of the spatial distribution of 
the particle emitting source because it is done implic- 
itly by the ellipse parameterization. A profile, linear in 
transverse rapidity is used in the filled ellipse case: 

p(r, 06) = (po + P2 cos{2(j)b) + p4 cos(4(/)b)) r{(j)b). (A5) 



The flow Fourier coefficients are defined by 

w„ = (cos [n((?!)p - «-)]), (A6) 

where (f>p is the azimuthal angle of the particle momen- 
tum. Assuming a Boltzmann plus flow distribution and 
longitudinal boost invariance, leads to the following ex- 
pression for w„: 



^ d<P, rdr /; d0pi^i(/3(r, q^b)) cos(n0p)e"(--^^) -o-i^^-'P.)n{r, <^,) 
"""^^'^^ /:^d0,/,°°rdr/;d0,Xi(/3(r,0,))e"('-.^^)-('A'-^.)fi(r,0,) ' ^ ' 

with a{r,(j)b) = {pt /T) smh{p{r, (j)b)) and j3t{r,(j)b) = (jnt/T) cos]i{p(r, (j)b)). The relation between (j)b and (j)s is given 
by Eq. (A3). All the integrals are done numerically in the filled ellipse calculation in order to preserve the possibility 
of computing interferometry radii, even though the formula can be simplified to: 

J^^d(f)b rdrKi{P{r,(j)b))cos{n(j)b)In{a{r,(j)b))n{r,(pb) 
"""^^^^ " J% d^b /o°° rdrK^r, cl,b))Io{a{r, MMr, d,b) ' ^ ' 



For the surface shell case the integral over r is trivial. 

APPENDIX B: DERIVATION OF THE MIXED 
HARMONIC EVENT PLANE METHOD 

^;i{EPi,EP2} 

Following the discussion in [23], we try to reduce the 
nonflow contribution of the first harmonic signal, vi, by 



subtracting the contributions to the flow vector perpen- 
dicular to the reaction plane from the component within 
the reaction plane. As an estimate of the reaction plane 
we use the second order event plane \E'2. Correlating the 
azimuthal angle of a particle, 0, with the first order event 
plane, ^'i, one then obtains 



(cos(0 - 'I'2) • cos(*i - *2) - sin((/) - ^'2) • sin(vl'i - '3>2)) = (cos((^ -I- v&i - 2*2)) = 
(C0S((/) - *Rp) • C0S(*1 - *Rp) • cos[2(*2 - *Rp)]) = 
(cos((/) - XE-Rp)) • (cos(*i - ^Rp)) • (cos[2(*2 - *Rp)]) = vi ■ Res(*i) • Res(*2)- 

I 



(Bl) 
(B2) 



The factorization in Eq. (B2) left-hand side is valid due to the resolution of the second order event plane, Rcs(vE'2)5 
the statistical independence of the three factors. While 
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can be obtained by calculating the square-root of the 
correlation of two subevent planes, the resolution of the 
first order event plane, Res(^'i), can be calculated by 
considering 



(cos^(*i 



(cos[2(^'i - *2)]) 

*R,p) • C0S[2(^'2 - *Rp)]) 



(C0S(*1 - *Rp))2 . (C0S[2(*2 - *RP)]) = 

Rcs2(«'i) •Res(*2)- 



(B3) 



Combining Eqs. (B2) right-hand side and (B3) yields 



z;i{EPi,EP2} = 



(cos((/. + ^1 - 2^2)) 
v/(cos[2(*i -*2)]) •Rcs(*2)' 



This approach is similar to the three-particle correlation 
method of Borghini, Dinh, and OUitrault [27]. One ob- 
tains their result by replacing the event plane angles ^"1 
and ^'2 in the right-hand side of Eq. (Bl) by emission 
angles of two particles [23]. 

Experimentally one wants to optimize the resolution of 
the second order event plane by measuring it in a region 



c where the signal of V2 is strong. This will be around 
mid-rapidity, preferentially. On the other hand, the influ- 
ence of nonflow can be reduced even further by measuring 
the azimuthal angle of the particle in one subevent, (j)"", 
and correlating it to the first order event plane in the 
other subevent, These subevents might by chosen 
randomly, or by dividing the acceptance into different 
regions in pseudorapidity. Since only half of all particles 
are used to determine each 4'° and the statistical 
errors are increased by a factor of \/2 compared to the 
three-particle cumulant method ?;i{3}. The final observ- 
able looks like this: 



^;i{EPi,EP2} = 



(cos(</>'^ + - 2*1)) 



v/(cos(«'J + - 2*^)) • Res(^'^) 



In our case each particle azimuth was correlated to the 
first order event plane determined in the other subevent 
within the FTPCs, and to the second order event plane 
measured in the TPC. 
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